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A class of discrete flavor-symmetry-based models predicts constrained neutrino mass matrix 
schemes that lead to specific neutrino mass sum-rules (MSR) . One of these implies in a lower bound 
on the effective neutrinoless double beta mass parameter, even for normal hierarchy neutrinos. Here 



t, we propose a new model based on the S4, flavor symmetry that leads to the new neutrino mass 

^_^ sum-rule and discuss how to generate a nonzero value for the reactor angle 613 indicated by recent 

^H experiments, and the resulting correlation with the solar angle 612 ■ 

PACS numbers: ll.30.Hv, 12.15.Ff, 14.60.Pq 



I. INTRODUCTION 



The discovery of neutrino oscillations have provided a strong evidence of the non- vanishing neutrino masses, although 
their nature (if they are Dirac or Majorana particles) has so far remained elusive. The observation of neutrinoless 
'n double beta decay (Oi//3/3) would provide an irrefutable confirmation of the Majorana nature of neutrinos [Ij. Majorana 
f^ neutrinos are characterized by a symmetric mass matrix whose parameters are restricted by the experimental data: 
C^ the neutrino oscillation parameters, such as mixing angles and neutrino squared mass differences [51|31, as well as the 
fvq limits on the QvjSp effective mass parameter [3i[5]. Indeed, upcoming Oiy/3/3 experiments are expected to improve the 

. . sensitivity by up to about one order of magnitude [BHS] . 
. !_( The general neutrino mixing matrix containing the three mixing angles and the CP violating phases can be 

rS parametrized in different equivalent ways |10H12| . A particular ansatz of the mixing matrix is the Tribimaximal 

^ Mixing Matrix (TBM) 13J which, despite the fact of the non-zero value of the 6*13 angle indicated by recent exper- 
iments [m [15] , can still be used as a good first approximation. Specially so taking into account that it can receive 
corrections from charged lepton diagonalization and/or from rcnormalization effects, depending on its scale of validity. 
Several flavor models based in non-Abelian discrete symmetries predict a two-parameter neutrino mass matrix 
which imply a particular mixing matrix form, as is pointed out in [13^ In Ref. |16) it was noted that in these models 
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only the following mass relations can be obtained, 



JL + ^^l^ 
7712 m'^ m\ ' 



1 



(1) 
(2) 

(3) 
(4) 



where x and ^ are free parameters that characterize each specific modeL A classification of all models predicting 
TBM mixing which generate mass relations similar to the first three are given there. The last case is completely new 
and here we will present a model from first principles, implementing the inverse seesaw mechanism [17| I18j as well as 
a non-Abelian flavor symmetry [T^, along the lines of Ref. [201, but adopting 5*4, instead of A4. 

The inverse seesaw scheme constitutes the first example of a low-scale seesaw scheme [H] with naturally light 
neutrinos. The particle content is the same as that of the Standard Model (SM) except for the addition of a pair of 
two component gauge singlet leptons, Vi and Si, within each of the three generations, labeled by i. The isodoublet 
neutrinos f^ and the fermion singlets Si have the same lepton number, opposite with respect to that of the three 
singlets i/f associated to the "right-handed" neutrinos. In the t^, t^'^, S basis the 9x9 neutral lepton mass matrix M^ 
has the form: 
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(5) 



where Y^ and Af are arbitrary 3x3 complex matrices, while fi is symmetric due to the Pauli principle. Following the 
diagonalization seesaw method in [22) one obtains the effective light neutrino mass matrix as: 



mJ,M'^^^fiM^^mD, 



(6) 



with the entry /i being very small. The diagram illustrating the mass generation through the inverse seesaw mechanism 
is shown in Figure [11 
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Figure 1: Inverse seesaw mechanism. 

It is straightforward to show that if mo and fi are both proportional to the identity, and 

' X y y 1 

M ^ Mtbm ^yx + zy — z, (7) 

\y y- z x + z J 

in the basis where the charged lepton mass matrix is diagonal. Here there is a specific (complex) relation among the 
parameters x, y and z [2^, leaving only two free complex parameters, we obtain the sum mass rule in Eq. (H). 



In the next section we give our model, in section III we present the predictions regarding the lower bound on the 
neutrinoless double beta decay amplitude, and discuss possible departures from tribimaximality, including a finite 
6*13 value. In the appendix [C] we present details on the symmetry structure, Yukawa couplings as well as the scalar 
potential of the model. 

II. THE MODEL 

Here we follow table I given in reference pO;, where some possibles schemes realizing the tri-bimaximal mixing 
pattern are summarized for the inverse seesaw case. From these we will implement case 1), since the other two cases, 
2) and 3), correspond to the mass sum rule relations C) and A) respectively, which have already model realizations 
in the existing literature, 

f A 0\ 
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In contrast to Ref. [501 here we adopt the S'4 flavour symmetry, instead of A^. 

In order to obtain the S'4-based inverse seesaw model we assign the charge matter fields as in table |l] Three right 
handed neutrinos vr are introduced, and three SU(2) fermion singlets S*;, i — 1, 2, 3, the latter transforming as the 3i 
(note that v'^ and vb. are conjugates, hence have opposite lepton number). All fermion fields in table ^ transform as 
the triplet 3i and the Higgs doublet as the trivial singlet li. 

On the other hand, in order to generate the desired mass matrix structures we introduce five flavon fields, 4>^, 0'^,, 
(j)i, (/}[, (f>i supplemented by the extra symmetries Z^ and Z2 , whose assignments are given in table Wn The presence of 
these extra Abelian symmetries in the theory ensures the presence of adequate zeros in the neutrino mass matrix. We 
keep renormalizability of the Lagrangian by adding a Frogatt-Nielsen fermion x and its conjugate x^, both singlets 
under the weak SU(2) gauge group [2^27] . In table In] we present the relevant quantum numbers of the matter fields 
in the theory under these extra symmetries. 
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Table I: Fields and transformation properties under 517(2), the 5*4 flavor symmetry, and global lepton number Ui{l) 
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Table II: Fields and their transformation properties under the Z3, and Z2 flavor symmetries 



The renormalizable Lagrangian relevant for neutrinos is 



Cy = YDi^LiVR^h + YXjVRiSj(f>y^ + Y^^^VR-Sjcf)'^ + fi^jSiSja^ 



(9) 



while the renormaUzable Yukawa terms involving the messenger fields is 

£x = M^xx' + Lhx + x'h'^i + X'h'^'i + X'h'^i- (10) 

After integrating out the messenger fields x, the effective Lagrangian for charged leptons takes the form 

(11) 



Ci = ^iLlR)h<j)i + ^{Lh)h(b'i + ^J^{LIr)H'{ 



where A is the effective scale. This effective Lagrangian is responsible for charged lepton mass generation, as shown 
in figure [2] 



■>y^i^i:^i 



X 



'li 



Figure 2: Diagram illustrating charged lepton mass generation. 



In order to obtain the desired neutrino mixing matrix we require the flavon fields to have the following alignments: 

((/.,)=«,(!, 0,0), {(t>i)=vi{l,l,l), {^i,)= VI, {1,1,1), (12) 

where we also define {(p,^') = Wj/', {4>i") = w/", {<j) ~ Va and {h) = v. In app endix [C| we report the form of the potential 
( C2 ) . We have verified that there exists a large portion of the parameter space where the required alignment is found 



to be a solution of minimization of the potential. 

With these alignments the three 3x3 blocks in Eq. ^ and charged lepton matrices take the form 

f ^J.v^ \ f Ydv \ f Y^,v^, \ 
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The charged lepton mass matrix, Eq. ( 14 ) , is diagonalized by the " magic" matrix 

/l 1 1 \ 



[/,.,= 
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1 a; w^ 

\1 OJ^ UJ j 



(13) 



(14) 



(15) 



On the other hand, by using Eqs. ([5]) and (rol it is straightforward to obtain the hght neutrino mass matrix which 
takes the form 

/^ \ 
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(16) 



where a = Y^iv^' /{y/Jw^Yj^v) and b — Y^v^ / {y/Jw^Y£,v) . In the basis where charged lepton mass matrix is diagonal, 
the Hght neutrino mass matrix is diagonahzed by the TBM-form, and the corresponding eigenvalues are given by 

1 



TO3 



(a+6)2 ' 

1 
(a-fc)2 ' 

1 



With these eigenvalues we obtain the neutrino mass sum rule 

1 2 
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(17) 



(18) 



which is, indeed, of the type given in Eq. Bl. 



III. PHENOMENOLOGY 



A. Neutrinoless double beta decay 



Using the symmetric parametrization of the lepton mixing matrix [101 lllj we can obtain the general expression of 
the mass parameter \mee\ which determines the Oi'/?/? decay amplitude as 
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(19) 



(PDG) , 

(symmetrical) . 

where c^ = cos 9ij and Sij — sin 9ij, mi,i — 1,2,3 the neutrino masses, and we adopt the symmetric parametrization 
where (j)i2 and 0i3 are the two Majorana phases. 

By varying the neutrino oscillation parameters in their allowed range one can plot \mee\ in terms of the lightest 
neutrino mass. Depending on which is the lightest neutrino one can have two different spectra, normal and inverse 
hierarchy, respectively. In the latter case one has a lower bound, on quite general grounds, as in this case there can 
be no destructive interference between the light neutrinos. 

In the present scheme, however, as noted in Ref. [inj, the neutrino mass sum-rule can be interpreted geometrically 
as a triangle in the complex plane, its area providing a measure of the Majorana CP violation. Then, fixing the (^, x) 
parameters for each model one can, in principle, determine the two Majorana CP violating phases. 

As a result there is a lower bound on \mee\ even in the case of normal hierarchy (for other schemes of this type see, 
for example, Ref. [TS] and references therein), as illustrated in Fig. [s] 



B. Quark sector 



Quarks are introduced as in table |IIl| where, differently from the lepton sector, we assign the fist and second 
families to a doublet representation of 5*4 and the third family to a singlet of S4, namely Qd = {Qi,Q2) ^ 2, 
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Figure 3: |mee| as a function of the lightest neutrino mass corresponding to the mass sum-rule in Eq. Q. The grey and dark 
bands correspond to generic normal and inverse hierarchy regions, while the green and yellow bands correspond to our flavour 
prediction varying the values of oscillation parameters in their 3 a C.L. range. The thin red bands correspond to the TBM 
limit. For references to the experiments see [6H9ll28j 



QRo = {QRn<lR2) ^ '^T Qs ^ Qa ^ li and = qu,^ ^ li. We add flavons ipD.s in doublet and singlet representations of 



the 5*4, 
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Table III: Quark sector and their transformation properties under the Z3, and Z2 flavor symmetries 



< - (Yi^QsdRs^'s + Y^Qo<i>DdR^ + YiQj,dR^<i>s + Y^Qj^^odBs + Y,^Qs^DdR^)h/A + h.c, (20) 



-C^ = iYrQsURs4>s + Y^Qo'^DUR^ + Y^QoUrJs + Y^Qd^durs + Y^^Qs'^oUR^h/K + h.c, (21) 

As in the charged lepton sector the dimension five operators can be given in terms of renormalizable interaction by 
introducing suitable messenger fields. Taking the VEV of (po in the direction (we verified that it is a possible solution 
of the potential) 



'>D 



^3,1) 



(22) 



the mass matrix for the quarks is 
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(23) 



which is very similar to the one proposed in Ref. 
and we refer to that paper for more detail. 



?5] where a fit of the quarks masses and mixing has been performed 



Finite 6*13 value 



As we have discussed so far the model leads to the tribimaximal mixing pattern. However by coupling an extra 
S'4-doublet flavon field one can obtain corrections from the charged lepton sector which induce nonzero values of ^13 
as recently suggested by the T2K [14] and Double-Chooz [15] results [2] including also recent reactor flux calculations. 

For example, consider a flavon scalar doublet under 6*4, <1> '^ 2 and transforming as (w, +) under Zj, x Z2. In the 
Lagrangian we must then include the term 



{LlR)h^. 



(24) 



This is a dimension five operator which can be obtained from a renormalizable Lagrangian by means of the messenger 
fields XtX'^ of table IT] as shown in figure (pi). Assuming that $ acquires VEV (<&) = (ui,U2), a natural vacuum 



alignment is ui — — v3u2, since this is consistent with the previous alignments in Eq. (12). Using multiplication rules 



in Appendix A one finds that the contribution from this term to the charged lepton mass matrix is 

/ -xl\vu2 \ 



5Mi = 



\l\vU\ + \j\vU2 



V 







\vu\ 
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(25) 



which modifies the diagonal entries SMi in the charged lepton mass matrix, Mi, so that the total Mi + SMi is 
no longer diagonalized by Ui^. This way one can induce a potentially "large" value for ^13, as hinted by recent 
experiments [TH [TS], and also potential departures of the solar and atmospheric angles from their TBM values. 
Moreover, in the presence of a nonzero ^13 one finds correlations among these neutrino mixing angles. The most 
interesting of these involves the solar angle, as illustrated in Fig. |4] The horizontal dot-dashed line indicates the best 
fit in the global analysis in Ref. [2]. On the other hand the region in the vertical band delimited by the dotted lines 
corresponds to the 3cr region found in Ref. [2] . One sees how the departure of the solar 6112 from its TBM value can 
be substantial and Iarge6'i3 values require the solar angle to lie below the TBM prediction ^. A nonzero ^13 would 
also open the way also for the phenomenon of CP violation in neutrino oscillations, one of the central goals of the 
upcoming generation of long baseline oscillation studies [30l [31] . 



^ For technical simplicity we have varied the charged lepton masscss well above what is allowed by their current determinations. The 
real correlation is expressed as a subband of this 
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Figure 4: Correlations between reactor and solar neutrino mixing angles. See text for explanations. 
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Appendix A: 54 group 

The Si group is the discrete group given by the four objects permutations. It contains 24 elements and can be 
obtained from two generators, S and T, satisfying: 



S'^^T^^ 1, 



ST'S = T. 



(Al) 



The Si irreducible representations are two singlets, li,l2, one doublet, 2, and two triplets, 3i,32. The product 
rules are given by (for more details see [32j.) 



Ij X Ij — l(i+j) „,od2 +1 V ?, J 

2 X 1, = 2Vi 

3j X Ij = l(i+j) ,„od2 +1 V i,j 

3i X 2 = 3i + 32 V i 

3i X 32 = I2 + 2 + 3i +32 

2x2:=li+2 + l2 

3i X 3, = li + 2 + 3i+32Vi 



(A2) 



where we can introduce the notation [fi x fi] and {/^ x /i} for the symmetric and antisymmetric part oi fj,Xfi respectively: 



[2 X 2] = li + 2, 

[3, X 3,] = li + 2 



{2x2} = I2 

h3i, {3, x3j = 32Vi 



(A3) 



Given the following representations 
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As the representation matrices can be choosen all reals, the conjugate representations transform in the same way. 
For the product of one dimensional representations the Clebsch Gordan coefficients are the trivial products of 
representations, and also for the product of the li singlet with any non-trivial representation 
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32. 



For the product with the I2 singlet 
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The Clebsch Gordan coefficients for the product 2x2 are 
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For the couplings 2 x 3i and 2 x 82, we have respectively 



' 0261 \ 

- 1(^80162 + a2&2) 
\ |(\/3ai63 - 0263) / 
' aibi ^ 

5(^80262 - 0162) 
\ -5(\/3a263 + 0163) J 



( 



3i, 



fllCl 



i(\/3a2C2 -aiC2) 

\ -|(\/3a2C3 +aiC3) j 

a2Ci 

-i(\/3aiC2 +0202) 

y 5(\/3aiC3 -0203) y 



\ 



3i, 



3,. 



And finally, for the 3i x 82 product 
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(All) 



Appendix B: Yukawa couplings 



Each term in the Lagrangian for neutrinos and charged Icptons can be decomposed by components as: 

liS ■ Sa ~ ^J-{SiSi + S2S2 + S'3S'3)cr, 



(Bl) 



YdL ■ vjih = YoiLiiyiR + L2V2R + LzVs.R)h, 



(B2) 



Y^{vr ■ S)(f)^ = Yi,[{u2rS3R + V3rS2)cJ)u1 + {'^IrSsR + i^3fl'S'l)0,.2 + {virS2R + V2RSl)(f>,y3], 



(B3) 



Yu'{vR ■ S)(t>,y' = Yy'{viRSi + V2RS2 + V-iRS3)(t>u' , 



(B4) 



Vl - Vl - - - - - - 

— {LlR)h(j)i = — [(£2^3fl + L3hR)h4>ii + {L1I3R + L3liR)h(j)i2 + {L1I2R + L2^i_r)/i0;3] 



(B5) 



V - V - - 

-j-{LlR)h(t>i> = -^[(L3^2fl - L2hR)h(j)l>l + {L1I3R -~ L3liR)h(t)i>2 + {L2I1R - il^2fl)/i0/'3] 
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^{LlR)h(f>l„ = |-(Li/iK + L2I2R + UhR)h(j,in 
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Appendix C: Scalar potential 

The most general renormalizable scalar potential is (without write the ^4 products explicitly): 

V = V{h) + V{a) + V{^,) + V{^,:) + V{^i) + V{<^v) + V{^v,) (CI) 

+y($,, $,s $/, $is $r') + Via, h, $„, $,.,$,, $r, $i"), (C2) 

with 

V{h) = ^ihh)h + \h{h)h){h)h), 
V{(j) = ;U^crV + Ao-(crV)(crV) 

i 

i 
i i 

i i 

i 

Via, h, $,, $,,,$,, $p, $,.0 - \''-{h^h)(a^<j) + X-'^^l^,){a^a) + X"'- {^1^,,){<J^ a) 

+A''^($|<i>0(fTV) + A''-($t cO;,)(ata) + A'"'^($|„$,,0('^^^), 
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i i 

i i 

t -i i 

-i i 

i i 

i i 

+ Yl ^«'''''' " {*r*r*i'*i" } . + H ^i" "' " {*|*J*i"*i" } 
+ ^Af'''"'"|$3',$3',$r *("}. + XI ^r"" 1*1*1" *i*i"}. 

i i 

+ Y ^^'"' {*;*(*.* J. + 51 ^"'"" {*/*p*.*.}. 

+ Y ^^'''"'^ {*^'*^'*-*-}^ + E ^r^" {*i*i"*.*.}i 

+ J2 ^^''"'" {*/'*/"*.*4» + I] ^P""" {*("*;"*.*.}, 

+ J2 Af " {*j*i*i:*.}^ + 51 Af "" {*J*;'*t*.}^ 

+ Y ^P'^'^ {*|'*p*t*4^ + I] K'""" {*|*r *t*4^ 

i i 

+ X Af '^''' {^l^V<^u^A^ + E ^^'''''' {*''*;'*-*-'}. 

i i 

+ E Af '"'''^' {$p$p.$,$,4^ + E Af''' {$3$z$t $^,|^ 

i i 

+ Y^f""' {*|*F*t*.'}^ + 5]Af"'''' {<f>J,<f>,$t$,,}^ 

i i 

+ Y ^^''''"' {*l'*i'*t*-' }^ + II ^»'"""' {*|*p'*t*.' } . 

i i 

+J2^r'' {*t*.*t*.'}, + E^r'"^''' {*t*.'*t*.'}, 

i i 

+5]Ar''''''{$t$,$t,$,,|^ 
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i i 

i i 

i i 

+h.c., 



where ^^ Xi {}^, ^^ k^ {}j sums over all possible ways to group the fields inside the brackets and make the product 
of representations in order to obtain a singlet. 
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